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Abstract 

We introduce a stochastic model of two-dimensional Brownian vortices associated with 
the canonical ensemble. The point vortices evolve through their usual mutual advection 
but they experience in addition a random velocity and a systematic drift generated by the 
system as a whole. The statistical equilibrium state of this stochastic model is the Gibbs 
canonical distribution. We consider a single species system and a system made of two types 
of vortices with positive and negative circulations. At positive temperatures, like-sign 
vortices repel each other ("plasma" case) and at negative temperatures, like-sign vortices 
attract each other ("gravity" case). We derive the stochastic equation satisfied by the 
exact vorticity field and the Fokker-Planck equation satisfied by the A-body distribution 
function. We present the BBGKY-like hierarchy of equations satisfied by the reduced 
distribution functions and close the hierarchy by considering an expansion of the solutions 
in powers of 1/N, where A is the number of vortices, in a proper thermodynamic limit. 
For spatially inhomogeneous systems, we derive the kinetic equations satisfied by the 
smooth vorticity field in a mean field approximation valid for A — > +oo. For spatially 
homogeneous systems, we study the two-body correlation function, in a Debye-Hiickel 
approximation valid at the order 0(1/N). The results of this paper can also apply to 
other systems of random walkers with long-range interactions such as self-gravitating 
Brownian particles and bacterial populations experiencing chemotaxis. Furthermore, for 
positive temperatures, our study provides a kinetic derivation, from microscopic stochastic 
processes, of the Debye-Hiickel model of electrolytes. 

1 Introduction 



Systems with long-range interactions are numerous in nature and share a lot of peculiar prop- 
erties p]. For example, it was recognized long ago that, for systems with unshielded long-range 
interactions, the statistical ensembles may be inequivalent at the thermodynamic limit (2J [3] . 
Therefore, we must carefully define the system under consideration before applying a statistical 
description. If we consider an isolated system of particles in interaction, where the microscopic 
dynamics is governed by Hamiltonian equations, the proper statistical ensemble is the micro- 
canonical ensemble. Alternatively, if we consider a dissipative system, where the microscopic 
dynamics is governed by stochastic Langevin equations, the proper statistical ensemble is the 
canonical ensemble. A general discussion of the dynamics and thermodynamics of Hamiltonian 



1 



and Brownian systems with long range interactions has been given in pQ and in [U El El HE] (de- 
noted Papers I-V) . To any Hamiltonian system with long range interactions, we can associate a 
Brownian system with long-range interactions by adding a friction force and a stochastic force 
in the equations of motion. For example, the self-gravitating Brownian gas studied in [9] is the 
canonical version of a stellar system [TU] and the Brownian Mean Field (BMF) model studied 
in [TT] is the canonical version of the Hamiltonian Mean Field (HMF) model [T2]. In order to 
emphasize the similarities and the differences between Hamiltonian and Brownian systems, it 
is of great conceptual interest to study these systems in parallel. 

The purpose of the present paper is to extend this study to the case of point vortices in 
two-dimensional hydrodynamics. Basically, a gas of point vortices constitutes a Hamiltonian 
system described by the Kirchhoff-Hamilton equations of motion where the coordinates (x, y) 
of the point vortices are canonically conjugate [131 EH- A statistical mechanics of this system 
has been developed by several authors [151 US El HH1 EE EUJ HH [221 123] following the pioneering 
work of Onsager [23]. For this system, the proper thermodynamic limit is such that the number 
of point vortices N — > +oo while their circulation tends to zero as 7 ~ 1/N and the area of 
the system remains finite V ~ 1. Assuming ergodicity, the system should reach, for t — > +00, 
a statistical equilibrium state described by the mean-field Boltzmann distribution. In order to 
vindicate (or not) this result and determine the timescale of the relaxation, we must develop a 
kinetic theory of point vortices. A general kinetic equation describing the collisional relaxation 
of point vortices has been derived in [25] from the Liouville equation by using a projection 
operator formalism. In a more recent paper [2B], we have re-derived this kinetic equation 
from a BBGKY-like hierarchy by using an expansion of the solutions in powers of 1/N for 
N — > +00, and we have shown that this equation is valid at the order 0(1 /N). In the limit 
N —* +00, the correlations between point vortices can be neglected and the kinetic equation 
reduces to the 2D Euler equation. This is the counterpart of the Vlasov equation in plasma 
physics. On the other hand, by considering an axisymmetric evolution, we have obtained an 
explicit expression of the collisional current at the order 0(1 /N) taking into account two-body 
correlations [251 126] . This leads to a kinetic equation that is the counterpart of the Landau 
equation in plasma physics. Collective effects can be included by developing a quasilinear 
theory of the Klimontovich equation [281 129] leading to a kinetic equation analogous to the 
Lenard-Balescu equation in plasma physics. It was numerically shown in [30] that the kinetic 
equation valid at the order 0(1 /N) does not in general relax towards the Boltzmann distribution 
because of the absence of resonances. This implies that the collisional relaxation time is larger 
than NtD (where to is the dynamical time). However, its precise scaling with N remains to 
be determined. This demands to develop a kinetic theory at the order 1/N 2 or smaller by 
taking into account three-body or higher correlations. This extension has not yet been done. 
Therefore, for the moment, there is no theoretical proof (coming from the kinetic theory) that 
the Hamiltonian point vortex gas relaxes towards the Boltzmann distribution for long times. 

In this paper, we shall study a different model of point vortex dynamics. Our aim is to 
introduce a stochastic model of point vortices associated with the canonical ensemble. In other 
words, we wish to construct a dynamical model of point vortices that leads to the Gibbs canon- 
ical distribution at statistical equilibrium. These vortices will be called Brownian vortices to 
distinguish them from usual vortices described by Hamiltonian equations (leading to the mi- 
crocanonical distribution at statistical equilibrium). By analogy with usual Brownian theory, 
if we want to introduce a canonical description, the first idea is to modify the usual point 
vortex model by adding a stochastic term in the Kirchhoff-Hamilton equations of motion. This 
extension has already been considered by Marchioro & Pulvirenti [31] in a different context. 

1 A similar approach has been developed simultaneously in |27| . 
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Their motivation was to take into account viscous effects in the point vortex dynamics. They 
considered a system of point vortices evolving through their usual mutual advection and ex- 
periencing, in addition, a random velocity whose strength is proportional to the square root 
of the viscosity. They showed that, in a mean-field limit valid for N — > +00, the evolution of 
the smooth vorticity field is described by the Navier-Stokes equation (instead of the 2D Euler 
equation). This model introduces a source of dissipation but it does not lead to the canoni- 
cal distribution at statistical equilibrium (the Navier-Stokes equation leads to the trivial state 
uj = for t — > +00). If we want to obtain the canonical distribution, we need to introduce, 
in addition, a drift velocity produced by the vortices as a whole. This drift velocity, which is 
perpendicular to the advective velocity, is the counterpart of the dynamical friction in usual 
Brownian theory. The drift coefficient, similar to a mobility, is determined so as to yield the 
Boltzmann factor at statistical equilibrium. It is related to the viscosity and to the inverse 
temperature by a sort of Einstein relation. 

This model of Brownian vortices was introduced in a Proceedings paper [32] and is here 
systematically developed. In a sense, this stochastic model takes into account the coupling 
with a thermal bath imposing the temperature rather than the energy. Therefore, it describes 
a system of stochastically forced vortices in contrast to the usual point vortex model describing 
an isolated system. At the present time, it is not clear how such a thermal bath can be realized 
in nature. However, this stochastic model is well-defined mathematically and constitutes by 
definition the dynamical model of point vortices associated with the canonical ensemble. There- 
fore, its study is interesting in its own right since it provides the out-of-equilibrium version of 
the situation described by the Gibbs canonical measure at equilibrium. As we said previously, 
it is of conceptual interest to compare the microcanonical and the canonical dynamics and 
study them in parallel. On the other hand, by a proper reinterpretation of the parameters, 
this stochastic model is isomorphic to the overdamped dynamics of Brownian particles with 
long-range interactions such as self-gravitating Brownian particles [9], bacterial populations 
experiencing chemotaxis [33], Debye-Huckel electrolytes [31] or electric charges moving in a 
strong vertical magnetic field (see Appendix [SJ. Therefore, it provides an interesting model 
of random walkers in long-range interactions and enters in the general class of stochastic pro- 
cesses studied in [HE]. In addition, this model allows for many generalizations with respect to 
previous works: (i) Point vortices can have positive and negative circulations. This generalizes 
the studies of self-gravitating Brownian particles and bacterial populations where the mass of 
the particles is always positive, (ii) Since point vortices have no inertia, the temperature can be 
either positive or negative. At negative temperatures, point vortices of the same sign have the 
tendency to attract each other (like stars in a galaxy □) and at positive temperatures they have 
the tendency to repel each other (like electric charges in a plasma). Therefore, by changing 
the sign of the temperature, we can study either attractive or repulsive interactions, (iii) The 
equilibrium states of point vortices in the canonical ensemble have been studied in detail and 
present peculiar features associated with the existence of critical temperatures beyond which 
there is no equilibrium [361 ETJ [371 EH]- Our dynamical model will allow to explore this range 
of parameters where rich and interesting behaviors are expected to occur. Therefore, on the 
viewpoint of statistical mechanics, this stochastic model of Brownian vortices has very rich 
properties and deserves a specific attention. 

The paper is organized as follows. In Sec. [21 we consider a single species system of Brow- 
nian vortices. Starting from the Fokker-Planck equation satisfied by the iV-body distribution 
function, we derive the BBGKY-like hierarchy for the reduced distribution functions (see Sec. 

2 The numerous analogies between the statistical mechanics of 2D point vortices at negative temperatures 
and the statistical mechanics of stellar systems have been described in [55] , 
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12. 2p . We close this hierarchy by considering an expansion of the solutions in powers of 1/N 
for N ^> 1. For N — > +00 (see Sec. 12.31) . we derive the mean field Fokker-Planck equation 
satisfied by the one-body distribution function (smooth vorticity field). This is the counterpart 
of the Vlasov regime in plasma physics. For spatially homogeneous systems, we derive the 
equation satisfied by the two-body correlation function at the order 0(1 /N) and investigate 
the stability of the system (see Sec. 12.41) . This is the counterpart of the Debye-Huckel theory 
in plasma physics. In Sec. 12.51 we derive the stochastic equation satisfied by the exact vorticity 
field before averaging over the noise. If we average this equation over the noise, we recover 
the first equation of the BBGKY hierarchy. In Sees. 12.61 and I2.7[ we consider specific forms of 
potentials of interaction corresponding to the ordinary Newtonian (or Coulombian) potential 
in two dimensions and to the Yukawa (or Rossby) potential appropriate to geophysical flows. 
For spatially homogeneous systems, we determine the two-body correlation function and the 
caloric curve relating the average energy to the temperature. We specifically distinguish the 
case of positive and negative temperatures. In Sec. [3] we consider a multi-species gas of Brow- 
nian vortices and focus particularly on the two-species system where the vortices have positive 
and negative circulations. We derive the BBGKY-like hierarchy and consider the Vlasov limit 
N — * +00 for the one-body distribution functions and the Debye-Hiickel approximation (valid 
at the order 0(1/N)) for the two-body correlation functions. Finally, in Sec. HI we conclude 
and discuss some perspectives of our work. 



2 The single species system 
2.1 The stochastic equations 

The canonical distribution of a single species system of point vortices with individual circula- 
tions 7 is 

Pn = \e~"\ (1) 

where H = J2i<j l 2u ( r ii r j) + \l 2 Z)i=i u ( r *> r «) * s the Hamiltonian, Z(j3) = J e~^ H dxi...dr N 
the partition function and (3 = 1/T the inverse temperature which can be either positive or 
negative. The free energy is F(/3) = — (1/(3) In Z((3). We have used the general Green's function 
formulation of Lin [39] , which is valid in an arbitrary bounded domain %. For convenience, 
we shall note U = J2i<j u ( r ^ r j) + \ J2f=i v ( r i> r «)- We ask: can we formally construct a 
dynamical model of point vortices associated with the canonical ensemble? The answer is yes. 
The dynamical model of point vortices consistent with the Gibbs canonical distribution (CQ) at 
statistical equilibrium is defined by the N coupled stochastic equations 

_L = _ 7Z x ViU - ^ 2 V,U + V2vRi(t), (2) 

where % = 1,...,N label the point vortices in the system. These stochastic equations define 
our model of two-dimensional Brownian vortices. The first term, where z is a unit vector 
perpendicular to the plane of motion, describes the mutual advection of the point vortices. 
The last term is a stochastic velocity where Rj(t) is a white noise such that (Rj(£)) = and 
(R?(t)R?(t')) = 8ij8 a p8(t - t'), and v > is a diffusion coefficient which plays the role of 
a viscosity (see below). Finally, the middle term is a drift velocity produced by the system 
as a whole, where fi plays the role of a mobility in Brownian theory. The drift velocity is 

3 In [19], the Green function u(ri,Tj) is noted — G(rj,rj) and the function u(ri,rj) is noted — g(ri,Ti). 
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perpendicular to the advective velocity. This systematic drift is necessary to compensate the 
effect of the stochastic velocity and yield the canonical distribution at equilibrium. This is 
similar to the dynamical friction in Brownian theory. 

The model ([2]) regroups and generalizes several models of particles in interaction already 
introduced in the literature. For /j, = v — 0, we recover the usual Hamiltonian model [13] of 
point vortices @: 

dv ■ 

_^ = - 7 zxV i t/, (3) 

which is associated with the microcanonical ensemble [2U UHl [201 [26]. F° r A* = and v > 0, we 
recover the stochastic equations considered by Marchioro & Pulvirenti [3Tj : 



^2 = _ 7ZX ViU + VtoRiit). (4) 
at 

Finally, without the mutual advection (obtained formally by taking z = 0), the stochastic 
equations 

-± = -iij 2 ViU + v^Ri(t), (5) 

describe the dynamics of Brownian particles in interaction in an overdamped limit where 7 
plays the role of the mass m and v plays the role of the diffusion coefficient (in that case, 
the model is valid in d dimensions) [HE]. In the case of attractive (e.g., Newtonian) potentials, 
these stochastic equations describe, for example, a gas of self-gravitating Brownian particles 
[SI ED] or bacterial populations experiencing chemotaxis [JTJ H2J H3] . 

2.2 Ensembles average 

Using the same procedure as the one developed in Paper II, we can readily write down the 
Fokker-Planck equation associated with the stochastic process ([2]). The evolution of the N- 
body distribution function P/v( r i5 ...,r/v,£) is given by 

9Pn mX^v 9Pn 9 ( 9Pn 2 P 9u\ 

i=l i=l x ' 

where V» = — 7Z x dU/dri is the total advective velocity of point vortex %. For v = \i = 0, we 
recover the Liouville equation corresponding to a microcanonical description [26]. For v > 
and n 7^ 0, the stationary solution of the Fokker-Planck equation (Q is 



It cancels individually the r.h.s. (Fokker-Planck term) and the l.h.s. (advective term) of Eq. 
([6]). If we compare this expression with the Gibbs canonical distribution (CD), we find that the 
inverse temperature is related to the mobility and to the diffusion coefficient (viscosity) by 

= ^, (8) 
v 

which is the Einstein relation in the present context. Since the vortices have no inertia, the 
temperature can take positive or negative values. Considering Eq. ([2]) with the ordinary 



4 The Hamiltonian point vortex system is isomorphic to a two-dimensional guiding center plasma in which 
charged rods move across an external magnetic field B = Bz with the self-consistent E x B drift [15] • In 
this analogy, the circulation 7 plays the role of the charge e and the stream function ip the role of the electric 
potential <f> (see Appendix [K\. 
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potential of interaction (J55|) or (1701) . we see that at positive temperatures (implying \i > 0) 
like-sign vortices repel each other and at negative temperatures (implying /i < 0) like-sign 
vortices attract each other. Let us introduce the free energy F = E — TS" where E = (H) = 
f P N Hdri...dr N is the average energy and S = — f P/vlnP/ V (ir 1 ...<ir7v the entropy. We also 
consider the Massieu functional J = —(3F = S — (3E given by 

J[P N ] = - J P N \nP N dr 1 ...dr N - -(3 J P N Hdv 1 ...dr N . (9) 

We see that the Gibbs state ([I]) is a critical point of J with the normalization constraint: 
5J-aJ 5P N dT 1 ...dr N = 0. Furthemore, 5 2 J = -(1/2) f(5P N ) 2 /P N dr 1 ...dr N < 0. Therefore, 
if the Gibbs state exists (i.e. if the partition function exists), then it is the only maximum of 
the free energy functional J[Pn}- Furthermore, inserting the Gibbs state (TjQ) in the Massieu 
functional ([91), we recover the equilibrium value of the free energy J {(3) = In Z(f3). On the other 
hand, using the Fokker-Planck equation ([6]) we find that 



1 / dP N 2n dU\ . 
vP N V ovi dvij 



Therefore, if \i ^ and v > the Fokker-Planck equation ([6]) satisfies an ^/-theorem appropriate 
to the canonical ensemble: J > and J = iff P/v is the canonical distribution (Q. The free 
energy (Q is the Lyapunov functional of the Fokker-Planck equation ([6]). From Lyapunov's 
direct method, we conclude that if J is bounded from above (or, equivalently, if the partition 
function exists), the Fokker-Planck equation ([2D will relax, for t — > +oo, towards the Gibbs 
distribution (TjQ). 

Proceeding as in Paper II, we can easily obtain the exact hierarchy of equations satisfied 
by the reduced distributions Pj(r±, ...,Tj,t). To that purpose, we first note that the stochastic 
equation ([2]) takes the form of Eq. (11-146) if we make the substitutions m — > 7, -D* — > v and 

duij d'uij 1 du i:j du i: j 
drj dr^ ^7 dr^ dr^ 

where d'uij/ dr-j is a convenient notation. We also need to take into account the term | J2f=i v ( r i> ] 
in the Hamiltonian. With these substitutions, we can immediately transpose the results of Pa- 
per II to the present context. The equations for the reduced distributions of order j = 1, N 
are given by 



dPj d 



(12) 

and they form the BBGKY-like hierarchy associated with the Fokker-Planck equation ([6]). The 
equilibrium distributions are obtained by cancelling the terms in brackets and by replacing & by 
d (since the steady state cancels the advective term). They satisfy the equilibrium BBGKY-like 
hierarchy 



5 Since the temperature can be negative in two-dimensional point vortex dynamics, it is often more convenient 
to use the Massieu functional (which is the direct Legendre transform of the entropy) rather than the usual free 
energy. By an abuse of language, we shall often refer to J as the "free energy" . 
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They can also be obtained directly from the Gibbs canonical distribution ([TJ by writing 



dP; 



N 



Or i 



■tfPsg-. 

or i 



(14) 



and integrating over r,- + i...rjy. The exact first two equations of the BBGKY-like hierarchy are 



dPx 
~dt 



_d_ 
dri 



dPi 

I 

9ri 



/i 7 2 (iV 



— — P 2 dr 2 + -//7 P 1 —v(t 1 ,t 1j 
or i 2 OTi 



(15) 



dP 2 

at 



_d_ 
dri 



dP 2 o „ d'u 12 ,„ T . 9 



OTi 



oti 



— — P 3 dr 3 + -m P 2 —v(r 1 , n, 

OTi 2 OTi 



+ (!<-> 2). (16) 



The next step is to decompose the reduced distributions in the form of a Mayer expansion by 
introducing the cumulants (see Eqs (1-14) and (1-15) of Paper I). Inserting the decomposition 
(1-14) in Eq. (ilHjl . we obtain 



dP 1 
~dt ~ 



_d_ 
dr x 

-1) 



or i 

d'u 12 



1) 



d'u 



12 



Oti 



P 1 {r 1 )P 1 (r 2 )dr 2 



1 & 
P 2 ( r i, r 2 )dr 2 + -irfp x —v{j: x , n) 

OTi 2 OTi 



(17) 



where P 2 is the two-body correlation function. Next, inserting the decomposition (1-14), (1-15) 
in Eq. ( |T6|) and using Eq. ( 1TT1) to simplify some terms, we get 



dP[ 
dt 



_d_ 

OTi 



OTi 



/x7 2 ^i(ri)Pi(r 2 



12 



OTi 



/X7 P 2 ( r l> r 2 



12 



OTi 



d'u 13 

OTi 

d'u 13 



OTi 

+(iy-2)/i7 2 



P 1 (n)P 1 (r 2 )Pi(r 3 )rfr 3 + (iV - 2)/i 7 2 
P^(r 1 ,r 3 )Pi(r 2 )rfr 3 + (iV-2) / i7 2 



^-^P^(r 1; r 2 , r 3 )c/r 3 + -/ry^-^-^n, n 



d'u 13 

OTi 
<9'mi 3 
OTi 

OTi 



P 2 '(r 1 ,r 2 )Pi(r 3 )dr 3 
P^(r 2 ,r 3 )Pi(n)rfr 3 
h (1 2), 



where P 3 is the three-body correlation function. These equations are exact for all N, but the 
hierarchy is not closed. 

We now consider the thermodynamic limit iV — > +oo in such a way that the normalized 
inverse temperature rj = f3N r j 2 remains of order unity. This corresponds to a regime of weak 
coupling since the "plasma parameter" firfu ~ 1/N — > 0. It is convenient to renormalize the 
parameters in such a way that the individual circulations scale like 7 ~ 1/N (so that the total 
circulation T ~ N'-f is of order unity), the inverse temperature scales like (3 ~ N and the 
domain area scales like V ~ 1. We also assume that v ~ 1 so that the diffusive timescale is of 
order unity. The dynamical time is also of order unity since to ~ ~ V/Y ~ 1. With these 
scalings, we have Pi ~ 1, |r| ~ 1, u ~ 1 and \i ~ /3 ~ N. Now, by considering the scaling of 
the terms appearing in the different equations of the BBGKY-like hierarchy, we find that the 
cumulants scale like Pj ~ 1/N^~ X . In particular, P' 2 ~ 1/N and P 3 ~ 1/N 2 . Therefore, up to 
the order 1/N, the first two equations of the BBGKY-like hierarchy take the form 



dPi 
~dt 



_d_ 

OTi 

-^ 2 N 



^ + ^ 2 (N-1) 
or 1 

d'u 12 



d'u 12 

OTi 



Pi(ri)Pi(r 2 )rfr 2 



1 & 
P 2 '(ri, r 2 )dr 2 + -fij 2 P 1 —v(r 1 , n) 
or 1 2 or 1 



(19) 
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dPi d 



dt dr 1 



^P^P^) — 



+N^ 2 J ^^(r a) r s )Pi(n)drs 



+ (!<-> 2). (20) 



The hierarchy is now closed because the three-body correlation function can be neglected at 
the order 0(1/N). 

On the other hand, the average energy 



N 

E = (H) = ^ / l 2 u ij P N dr 1 ...dr N + - ^ / l 2 v ii P N dr 1 ...dr N , 

i<j J i=l J 



(21) 



can be written 



1 f NT 

E = -N(N - 1) 7 2 / P 2 (ri, r 2 )« 12 dr 1 dr 2 + y 7 2 / PifaMn, n)rf ri . (22) 

Substituting the Mayer decomposition (1-14) in the previous equation, and introducing the 
smooth vorticity field u = N^Pi and the stream function 



^(r,t) = J u(r,v')uj(r',t)dr', 



(23) 



we obtain 



E = \J ui/idr-N^ J P 1 (r 1 )u 12 P 1 (r 2 )dr 1 dr 2 

1 f N f 

+-N{N - 1) 7 2 / P 2 '(n, r 2 )u 12 dr 1 dr 2 + - 7 2 / Pi^Mn, n)dn. (24) 

2.3 The mean field approximation 

For N —* +oo, we can neglect the two-body correlation function P' 2 ~ 1/N. In that case, the 
two-body distribution function factorizes in two one-body distribution functions: 

P 2 (r l7 r 2 , t) = P 1 (r 1 ,t)P 1 (r 2 , t) + 0(1/ N). (25) 

More generally, for N — > +oo, we have the factorization 

AT 

Pv(ri, ...,rjv, t) = JJPi( ri ,t). (26) 

i=l 

This is the equivalent of the mean field approximation (valid in the Vlasov regime) in plasma 
physics or stellar dynamics. The mean field approximation is exact at the thermodynamic limit 
N — > +oo. Introducing the smooth vorticity field u = N^Pi and the stream function f[2"3"j) . 
the first equation fTT^l) of the BBGKY-like hierarchy reduces to the mean field Fokker-Planck 
equation |^ 

— + m-Vuj = uV- (Vuj + /3 7 cjV^) = -V • J, (27) 



6 To avoid confusions, it should be emphasized that the mean field Fokker-Planck equation l|27[) has a com- 
pletely different origin and interpretation from the Fokker-Planck equation (140) derived in [26] . In |26j . the 



s 



where u = — z x Vi/j is the mean field velocity. The mean field Fokker-Planck equation (1271) 
conserves the circulation T = N'y = J u)dr, proportional to the vortex number. On the other 
hand, its steady solution is the mean field Boltzmann distribution 



to r 



(28) 



This distribution cancels individually the advective term (l.h.s.) and the Fokker-Planck term 
(r.h.s.) of Eq. fT2~Tj) . The mean field Boltzmann distribution can also be obtained from the 
equilibrium BBGKY-like hierarchy in the limit N — > +00 (see Paper I). Substituting Eq. 
(I28p in Eq. (l23l . we obtain an integro-differential equation Vina; = — /?7 J Vu(r, r')u>(r') dr' 
determining the vorticity field at statistical equilibrium. 

For N — > +00, using Eq. f[2T)j) . the Boltzmann entropy S = — J Pjy In P N dri...dr N reduces 

to 



S = -N Pi(r,t) lnPi(r,t)dr. 



(29) 



It can be written (up to an additive constant) 



su 



— In — dr. 

7 7 



(30) 



This expression of the Boltzmann entropy can also be obtained from a classical combinatorial 
analysis [T5l l30] . On the other hand, for N —* +00, we obtain from Eq. ( l24|) the mean field 
energy 



EM 



- I uipdr. 



(31) 



The mean field free energy (more precisely the mean field Massieu function) is J[u] = S[uj] — 
0E[u] where S and E are given by Eqs. (|30|) and (13TT) . It can be obtained from the free energy 
([9]) by using the mean field approximation (|26|) valid for N — * +00 [5j. The kinetic equation 
( I27p can then be rewritten in the form 



du(r, t) 
dt 



z x V 



SE[u] 
5lj(t, t) 



Vw(r,t) 



SJ[u] 
Suj(r, t) 



(32) 



The steady state corresponds to a uniform value of a = 5J/8u(r) leading to the Boltzmann 
distribution (12811 . On the other hand, it is straightforward to establish that 



J 



-dr 



V^Vdr>0. 

ouj J 



(33) 



Therefore, if \i ^ and v > the mean field Fokker-Planck equation (1271) satisfies an H- 
theorem appropriate to the canonical ensemble: J > and J = iff u is the mean field 



Fokker-Planck equation (140) describes the relaxation of a test vortex (tagged particle) in a fixed distribution 
of field vortices at statistical equilibrium (thus tp(r) is a given function independent on time). Here, the mean 
field Fokker-Planck equation (|27|) describes the evolution of a Brownian system of point vortices as a whole for 
an assumed microscopic dynamics of the form ([2|) (thus ip{ r ,t) is a function of time produced self-consistently 
by the distribution of point vortices according to Eq. ([25)) ). The mean field Fokker-Planck equation l|27p is also 
physically different from the relaxation equation (12) with a time dependent temperature (3(t) introduced by 
Robert & Sommeria [44] to phenomenologically describe the violent relaxation of the 2D Euler equation on a 
coarse-grained scale. 
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Boltzmann distribution (1251) . The free energy J[u] is the Lyapunov functional of the Fokker- 
Planck equation (I2~T1) . The Boltzmann distribution (I28|) is a critical point of J at fixed T 
(cancelling the first variations) and it is linearly dynamically stable iff it is a (local) maximum 
of J at fixed T. If J is bounded from above, we conclude from Lyapunov's direct method, that 
the mean field Fokker-Planck equation (12 7p will reach, for t — > +00, a (local) maximum of J 
at fixed circulation. If several local maxima exist, the selection of the maximum will depend 
on a complicated notion of basin of attraction. If there is no maximum of free energy at fixed 
circulation, the system will have a peculiar behavior and will generate singularities associated 
with vortex collapse (see Sees. 12.61 and l2TT|) . We note that the relaxation equation (1271) can be 
obtained from a maximum entropy production principle (MEPP) by maximizing the production 
of free energy J at fixed circulation (and other physical constraints) [l5j 06]. 

Finally, we note that the corresponding equation for the velocity field u = — z x is 

+ (u • V)u = — Vp + uAu - vpjuu, (34) 
ot p 

where p(r, t) is a pressure field. Indeed, taking the curl of this equation and using ujz = V x u, 
we recover Eq. (l2~Tj) [the intermediate steps of the calculation make use of the identity of vector 
analysis Au = V(V ■ u) — V x (V x u) which reduces to Au = —V x (ujz) = z x Vcu for an 
incompressible velocity field u = — z x Vip. On the other hand, V x (z x Vcu) = Auz and 
V x (uu) = x u + Vo; x u = —V ■ (uVip)z; finally, the transformation of the advective 
term is standard]. Interestingly, we note that the systematic drift of the point vortices —fi'-fVip 
in Eq. (TSTj) corresponds to a friction force — fi'yuu in the equation (13~4"|) for the velocity. This 
"duality" is a particularity of the Brownian point vortex system which has no counterpart in 
systems of material particles. 



2.4 The two-body correlation function 



The evolution of the two-body correlation function at the order 1/N is determined by the second 
equation (120]) of the BBGKY hierarchy. Let us consider a spatially homogeneous distribution of 
point vortices (to leading order) so that Pi(r,t) = P + P(r,t) where P is of order 1/N. Here, 
we consider an infinite domain so that u(r, r') = u(\r — r'\). Writing P 2 '(r, r', t) = P 2 /i(|r — r'|, t), 
the second equation (120)) of the BBGKY hierarchy equation simplifies into 



dh 



2uA 



h + (3-i 2 u + (3n-i 2 / /i(y,t)u(|x-y|)dy 



(35) 



where x = r — r'. The steady distribution of this equation satisfies an equation of the form 



dh 
<9x 



du 
<9x 



j h(y,t)^(\x-y\)dy. 



(36) 



Noting that the integral is a convolution, this integro-differential equation is easily solved in 
Fourier space leading to 



1 + (2ir) 2 (3n-f 2 u(k)' 



(37) 



The dynamical equation (1351) can be solved similarly and the temporal evolution of the two-body 
correlation function in Fourier space with initial condition /i(k, 0) = is found to be 



h(k,t) = h eq (k)\l-e 



-2uk 2 \l+(2ir) 2 l3nr 



! «(fc)]t| _ 



(3* 
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From this equation, we see that the homogeneous equilibrium distribution is stable iff 



1 + {2-K) 2 (3n 1 2 u{k) > 0, (39) 

for any wave vector k (otherwise, the wave vectors that do not satisfy this relation correspond to 
unstable modes that grow exponentially rapidly with time). This inequality usually determines 
a critical temperature (3q (similar to a spinodal point) beyond which the homogeneous system 
becomes unstable [I]. The stability criterion f[3"9"j) can also be obtained by studying the linear 
dynamical stability of a spatially homogeneous steady solution of the mean field Fokker-Planck 
equation (j27j) with uj = N'yPo and u = 0. Indeed, decomposing the perturbation in normal 
modes as Scu ~ e l ' kT "' 7t ' 1 the dispersion relation obtained by linearizing the Fokker-Planck 
equation (T2T|) is ia = vk 2 [1 + (2w) 2 f3wy 2 u(k)}. The system is stable iff io > for all k 
corresponding to the inequality (I3"9"j) . This inequality is also implied by Eq. (11491) of Appendix 
|B] and by the condition that the equilibrium state must be a maximum of J[u] at fixed T[u] 
(see Sec. 4.4. of @]). 

2.5 The exact vorticity field 

In this section, we shall derive the stochastic equation satisfied by the exact vorticity field before 
averaging over the noise. To simplify the expressions, we shall work in an infinite domain where 
Uij = u(\ri — rj\) and Vu = 0. The generalization of our results to a bounded domain is obtained 
by replacing u(\ri — | ) by w(rj, i\,-) and by adding a term (l/2)7 2 Vw(r, r) in the velocity field. 
The exact vorticity field is expressed as a sum of Dirac distributions in the form 

N 

w d (r,t)= 7 X>(r -*(*)). (40) 

i 

Adapting the calculations of [HI HH HHl EHl El H3] to the present context, we find that the exact 
vorticity field is solution of the stochastic equation 



dt 



[r,t)-zx V- (co d {r,t)V J u d {r',t)u{\r-r'\)dr') = uAu d (v,t) 



+/i 7 V- U d {r,t)V J ^(r^Mlr-r'Odr'J + V - {^2u 1 uj d {v, t)R(r, t) 



(41) 



where R(r, t) is a Gaussian random field such that (R(r,t)) = and (R a (r, t)i? /3 (r / , t')) = 
<W£( r — r')<5(t — t'). Introducing the exact stream function 

M*,t) = J u d (r',t)u(\r-r'\)dr', (42) 

and the exact velocity field = — z x Vip d , the stochastic equation (T4T!) can be rewritten in 
the form 

^ + u d ■ Vuo d = uAu d + V ■ (ji-yuj d Vil; d ) + V • ( y/2u^u d R) . (43) 

Let us consider some particular cases. For v = fi = 0, we find that the exact vorticity field is 
solution of the 2D Euler equation 

^ + u d ■ Vuo d = 0. (44) 
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This is the counterpart of the Klimontovich equation in plasma physics. On the other hand, 
for the model considered by Marchioro & Pulvirenti [31] where \x = and v > 0, the exact 
vorticity field is solution of the stochastic equation 

duj d 



dt 



+ u d ■ Vu d = uAu d + V ■ (y/2h"yu d R). 



(45) 



This equation was not given in [31]. Introducing the exact free energy functional 

J = S-(?H = - [ ^^ln^^dr-i/3 J u d {r,t)u{\r-r'\)uj d {r\t)drdr\ (46) 



J 7 7 
the stochastic equation ( l43|) can be rewritten in the form 



-V 



du d (r,t) 
dt 
SJ[oJ d ] 



z x V 
V 



6H[u d ) 



y/2wyu d (r,t)R(T,t) 



(47) 



5u d (r,t)_ 

It can be shown [8] that this stochastic equation reproduces the equilibrium two-body corre- 
lation function (Su>(t)5u(t')) that is related to the function h(r — r') determined in Sec. 12.41 
according to the identity (11471) . Equation (JiTjl can be viewed as a Langevin equation for the 
field u) d (r, t). Using standard methods [50], the corresponding Fokker-Planck equation for the 
probability W[o;d, t] of the field u d is 

dw[oj d ,t] , f 5 |Y_ wr7 <wm 



3t 



5w d (r,t) 



z x V 



— 1/7 



5 



<W(r,t) 



V-w d (r,t)V 



<M r >*) 
5 tfJM 



V^(r,t)W^,t] jdr 



(48) 



_5^ d (r,t) 5w d (r,t)_ 

and the equilibrium distribution is ^[cu^] oc e J ^ d ^~ auJd . 

If we now average Eq. (14T1) over the noise, we find that the evolution of the smooth vorticity 
field u;(r, t) = (u d ) is governed by an equation of the form 

|(r,i)-zxV' / Mr, t)u d (r', t))Vu(|r - r'\)dr' 



uAu(r, t) + /17V ■ J (u d {r, t)u d (r', *))Vu(|r - r'\)dr'. 



(49) 



This equation is equivalent to Eq. ([15]) giving the exact evolution of the one-body distribution 
function. Indeed, using u(r,t) = N^fPi(r,t) and the identity (see Appendix IB"]) : 

Mr, t)uj d {v\ t)) = iV 7 2 Pi(r, t)5(r - r') + N(N - l) 7 2 P 2 (r, r', t), (50) 

we find that Eqs. ( 149]) and ( Tl~5]) coincide. Furthermore, if we make the mean field approximation 
(^(r, t)u>d( r ', t)) = u(r, t)ou(r', t) valid for iV — > +00, we obtain 

duo 



dt 



+ u ■ = V ■ (z/Vu; + /rywV^), 



(51) 



which is equivalent to Eq. (127j) . In Fourier space, the mean field Fokker-Planck equation (|5l 
can be written 

<9u> 



(k, *) + (2vr) 2 J k • k' ± o>(k - k', t)u(k')u>{k' , t)dk' 
-uk 2 u(k, t) - (2tt) 2 /i7 J k ■ k'cj(k - k', t)n(kXk', t)dk', 



(52) 
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where we have noted ki = z x k. Let us consider some particular cases of Eq. (|51|) . For 
v = \i — 0, we find that the smooth vorticity field is solution of the 2D Euler equation 

^ + u • Vlu = 0. (53) 

This is the counterpart of the Vlasov equation in plasma physics. On the other hand, for the 
model considered by Marchioro & Pulvirenti [31] where /i = and v > 0, the smooth vorticity 
field is solution of the Navier-Stokes equation 

^ + u ■ Voo = uAuj. (54) 
at 



The derivation of this mean field equation was proven rigorously in [3 1 J . 

Equation floTjl for the ensemble averaged vorticity field w(r,t) is a deterministic equation 
since we have averaged over the noise. In contrast, Eq. (J4TJ) for the exact vorticity field t) 
is a stochastic equation taking into account fluctuations. However, it is not very useful in 
practice since the field Ud{r, t) is a sum of Dirac distributions, not a regular function. Therefore, 
it is easier to directly solve the stochastic equations ([2]) rather than the equivalent Eq. (14T|) . 
Following [501 E] , we can keep track of fluctuations while avoiding the problem of (^-functions by 
defining a "coarse-grained" vorticity field uJ(r, t) obtained by averaging the exact vorticity field 
on a spatio-temporal window of finite resolution. For a weak long-range potential of interaction 
and for a sufficiently small spatio-temporal window, we propose to make the approximation 
uj( 2 \r, r', t) ~ u;(r, t)Zu(r', t). Then, we find that the "coarse-grained" vorticity field satisfies a 
stochastic equation of the form 



, u- VuJ = uAU + V ■ (fjrftjV^) + V • ( v / 2^tEJR). (55) 

If we ignore the noise term, Eq. (1551) reduces to Eq. (j27j) . In that case, the system tends 
to a steady state that is a maximum (global or local) of the free energy functional J[uJ] at 
fixed circulation (minima or saddle points of free energy are linearly dynamically unstable). 
If the free energy admits several local maxima (metastable states), the selection of the steady 
state will depend on a notion of basin of attraction. Without noise, the system remains on a 
maximum of free energy forever. Now, in the presence of noise, the fluctuations can induce 
dynamical phase transitions from one maximum to the other. Thus, accounting correctly for 
fluctuations is very important when there exists metastable states (see discussion in Sec. 2.3 
of [8]). This would be an interesting effect to study in more detail. 

2.6 Newtonian or Coulombian potential 

Let us make the previous results more explicit by considering particular forms of potentials of 
interaction. The usual potential of interaction between point vortices is solution of the Poisson 
equation 

Au = -6{x). (56) 

It corresponds to a Newtonian (or Coulombian) potential in two dimensions. For box-confined 
vortices interacting through this potential, it can be shown that the Gibbs canonical distribution 
exists (i.e. the partition function is finite) if, and only, if [211 E3 EH] 

P > A ^ (57) 
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At positive temperatures (3 > 0, point vortices of the same sign tend to "repel" each other and 
accumulate on the boundary of the domain. This corresponds to a repulsive interaction like 
between electric charges of the same sign in plasma physics. At negative temperatures f3 < 0, 
point vortices of the same sign tend to "attract" each other and form clusters. This corresponds 
to an attractive interaction like between stars in a galaxy. As the inverse temperature is reduced 
the cluster is more and more condensed and for (3 = /3 C , the equilibrium distribution is a Dirac 
peak containing all the vortices. For (3 < j3 c , there is no equilibrium state (the Gibbs canonical 
distribution (TjQ) is not normalizable and the free energy (Q has no maximum). This regime can 
be studied dynamically with the stochastic model (j2J) introduced here. 

In the mean field approximation valid for N — > +oo, the evolution of the smooth vorticity 
is described by the coupled system 

^ + u-Vlo = uV ■ (Voo + (3-fujVifj), (58) 



- Ai[) = co. (59) 

The equilibrium state is obtained by solving the Boltzmann-Poisson equation [35] . By a proper 
reinterpretation of the parameters, these equations are isomorphic to the 2D one-component 
Smoluchowski-Poisson (SP) system, describing a gas of self-gravitating Brownian particles [9], 
with an additional advective term. Like in gravity, we expect that for (3 < /3 C , this system 
will describe the collapse of the Brownian vortex gas. In fact, if we consider an axisymmetric 
evolution, the advective term cancels out in Eq. fl58l) and we can immediately transpose the 
results obtained in gravity for the SP system: (i) For (3 > f3 c , the system fl58l) - fl59l tends to an 
equilibrium state in a bounded domain [52] and evaporates in an infinite domain [53]. (ii) For 
(3 = /3 C , the system forms, in an infinite time, a Dirac peak containing the N point vortices. 
The central vorticity diverges exponentially rapidly with time in a bounded domain [52] and 
logarithmically in an unbounded domain [53]. (iii) For j3 < (3 C , the system forms, in a finite 
time t co u, a Dirac peak containing {(3 C / (3)N point vortices surrounded by a halo of vortices 
evolving pseudo self-similarly [52] . A Dirac peak containing all the point vortices is formed in 
the post-collapse regime t > t coU in a finite time t end [54"j . 

Let us now determine the two-body correlation function of an infinite and homogeneous 
distribution of point vortices using Eq. fl36|) . In fact, for (3 ^ 0, a single species system of point 
vortices with potential of interaction fl56|) cannot be spatially homogeneous. However, our 
approach can be justified by the following arguments: (i) We can proceed as if the system were 
infinite and homogeneous, making the equivalent of the "Jeans swindle" in astrophysics [55] . 
This can give us some hints about the behaviour of the correlation functions for our system, 
(ii) We can add a neutralizing background, like in the Jellium model of plasma physics, so 
that spatially homogeneous equilibrium states of point vortices now exist, (iii) We can consider 
a two-species system and focus on one of the two species. It will be shown in Sec. [3] that 
the correlation functions of a two-components, globally neutral, homogeneous system of point 
vortices are the same as those derived here by making the "Jeans swindle" . This remark justifies 
a posteriori the validity of the following results for a neutral system of point vortices. 

In an infinite domain, the potential of interaction and its Fourier transform are given by 

u{x) =~ ln(x), (2ir) 2 u(k) = 1. (60) 

Taking the divergence of Eq. ( 136]) and using the Poisson equation ( l56i) . we find that the two- 
body correlation function satisfies an equation of the form 

Ah-pn^h = /3 7 2 5(x). (61) 
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The Fourier transform of the correlation function is 



According to the criterion (|39|) . the system is stable for f3 > and unstable for /3 < 0. At 
positive temperatures, a point vortex of given circulation has the tendency to be surrounded by 
point vortices of opposite circulation (in the two-species gas). This is similar to Debye shielding 
in plasma physics and this allows the existence of stable homogeneous states. At negative 
temperatures, an infinite and homogeneous distribution of point vortices has the tendency to 
collapse and form clusters. This is similar to the Jeans instability in self-gravitating systems. 
Let us consider these two cases consecutively. 

At positive temperatures, it is convenient to introduce the Debye wavenumber ko — (f3n , y 2 ) 1 ^ 2 . 
The two-body correlation function can then be expressed as 

-k 2 Ft<y 2 
n(27r) 2 h(k) = h(x) = -^K (k D x). (63) 

This corresponds to the two-dimensional Debye-Hiickel theory. The correlation length is equal 
to the Debye length 

Ld = v D = wiw (64) 

The energy of interaction (correlational energy) is given by 

E = InVV / or E = ^n 2 ^ 2 V(2ir) 2 J h(k)u(k)dk. (65) 

Evaluating the first integral in Eq. (j65p . which is convergent, we get 



ln(^) + 2 7i 



P = -l_ e -^- 27E , (66) 



where je = 0.577... is Euler's constant 0. On the other hand, using Eq. fl63l) . the spatial corre- 
lations of the vorticity in position space and its Fourier spectrum are given by (see Appendix 
©: 

(Mr)Mr')> = n^Ur - r') - ^K (k D \v - r'|) j, (67) 



{6M&u)= ', ■ ^ (k + kO- (68) 



The spectrum is depressed at small k (i.e. long-wave fluctuations are reduced) corresponding to 
Debye shielding in plasma physics. The velocity distribution W(V) and the spatial correlations 
(V(0) • V(r)) of the velocity field created by the point vortices are discussed in [56]. 

7 The second integral in Eq. |65|) is divergent. If we introduce a lower cut-off k m i n ~ V^ 1 ^ 2 taking into 
account the large but finite extent of the system, we get E — -~n 2 j 4 f3V/(4:Tr) J k °° (l/fc 2 )(l/(fc 2 + k 2 D ))kdk. 

This leads to the following relation E = -N-y 2 /(8n) In (l + fiN^/Sn) or (3 = 8-K/(N-/ 2 )(e- 87rE ^ N ^ - 1) 
between the average energy and the inverse temperature. For convenience, we have taken k max = (Sn/V) 1 ^ 2 
in order to recover the critical temperature (i c = —8ir/(N^ 2 ) for E — > +oo. However, we recall that the above 
results are only valid for (3 > 0. 
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At negative temperatures, we can introduce the equivalent the Jeans wavenumber kj = 
(—Pwy 2 ) 1 / 2 . The two-body correlation function can then be expressed as 

<^fKk) = ^fpr, h(x) = ^Y (kjx). (69) 

We note that h(x) = —f3^ 2 u(x) for x — > 0. The oscillatory behaviour of the correlation function 
h(x), or the fact that the denominator of h(k) is negative for k < kj, is related to the instability 
(collapse) of a homogeneous distribution of point vortices at negative temperatures (see criterion 
( 139]) ). The system is unstable to large wavelengths, above the Jeans length Lj ~ fcj 1 . When 
the system becomes spatially inhomogeneous, the dynamics of clustering can be studied with 
the mean field equations (1551) - (I5TH) . 



2.7 Yukawa or Rossby potential 

We now assume that the potential of interaction is solution of an equation of the form 

Au- klu= -£(x), (70) 

where k^ 1 is the characteristic range of the potential. This screened Poisson equation describes 
geophysical flows in the so-called Quasi-Geostrophic (QG) approximation. In that context, ko 
is the Rossby wavenumber. The screened Poisson equation ( J7D1 also arises in the Hasegawa- 
Mima model [57] which describes low-frequency drift waves in magnetized plasmas. In an 
infinite domain, the potential of interaction and its Fourier transform are given by 

u(x) = ^K (k x), (27r) 2 u(k) = (71) 

For a spatially homogeneous system, using Eqs. fl37|) and flTTT) . the Fourier transform of the 
correlation function is 

According to the stability criterion fl39l) . the system is stable iff 

P > — % = A)- (73) 

We note that, for a screened Coulombian potential, stable homogeneous distributions exist at 
negative temperatures, above a critical temperature /3 depending on the characteristic range 
kg 1 of the potential. In terms of this critical temperature, the correlation function f)72p can be 
rewritten 

= k > + t l{?l my "(x) = Mi - ew*] • (74) 



The correlation length is 



L = l(l-/?//5 )- 1/2 . (75) 
k 



Let us first consider the case of positive temperatures. Introducing the Debye wavenumber 
kf) = (Pwy 2 ) 1 / 2 , the correlation function and the correlation length can be written 

n(2n) 2 h(k) = - , 2 , L = -=L=. (76) 



k 2 + k 2 + k 2 D ' MTk\ 



D 
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At positive temperatures, the correlation length is smaller than the Rossby length: L < L 
where Lq = 1/ko- This corresponds to the formation of a screening cloud of opposite sign 
vortices around each vortex due to collective effects (like in plasma physics). For T — > + , i.e. 
(3 — > +00, the screening length coincides with the Debye length: L = — > 0. For T — > +00, 
i.e. (3 — > + , collective effects are negligible and we have h(x.) = — /?7 2 m(x) [1] so that the corre- 
lation length coincides with the Rossby length: L = L . At negative temperatures, introducing 
the Jeans wavenumber kj = (— /Jrry 2 ) 1 / 2 , the correlation function and the correlation length 
can be written 

At negative temperatures, the correlation length is larger than the Rossby length: L > L . 
This regime corresponds to the formation of a reinforcing cloud due to the clustering of point 
vortices of the same sign (like in gravity). For T — > — 00, i.e. (3 — > 0~, collective effects are 
negligible so that the correlation length coincides with the Rossby length: L = L . For T — > T , 
i.e. /9 — > /So, the correlation length diverges. For (3 < (3q, the system is unstable to wavenumbers 

k< k max = k Q {f3/^-l) 1 / 2 . (78) 

For (3 = Po, we have /c maa; = and for (3 — > — 00, we have k max ~ /cj — > +00. For a screened 
potential, since fc^ aa , = kj — k 2 ,, the lengthscale at which the system becomes unstable is larger 
than the Jeans length. 

For a screened potential, the energy of interaction f!65p is given by 

jji n^(3V [ + °° 1 1 , 7Q , 

4tt io P + A^ + ^(l-/V/3o) ' 1 j 

yielding 



^7 1 /, nm^ P 



8tvE 



E = -—!-hi(l-l3/i3o), j- = l-e^. (80) 

87T /^o 



On the other hand, using Eq. (1741) . the spatial correlations of the vorticity in position space 
and its Fourier spectrum are given by (see Appendix IB"]) : 

(Mr)Mr')) = ^7 2 {^(r - r') - ^K [k (l - /3//3 ) 1/2 |r - r'|] 



(&S k ^)= , ^ J(k + kO. (82) 

Interestingly, these results can also be obtained directly from the density of states (for large 
negative energies) fTf\ [T9] or from the partition function [51] by making the random phase ap- 
proximation. The complementary approach presented here, based on an equilibrium BBGKY- 
like hierarchy, is simpler because it avoids the use of complex analysis. In the present context of 
Brownian vortices (canonical ensemble), the control parameter is the temperature (3 while for 
usual vortices (microcanonical ensemble) [HIIIH], the control parameter is the energy E. Let us 
interpret the results (!HT]) - (!82|) in terms of the temperature making a parallel with the discussion 
given by [TF] in terms of the energy. When T = 00 (i.e. (3 = 0), the spectrum (1821) is flat as 
for a random distribution of vortices. When (3 > 0, the spectrum is depressed at small k (i.e. 
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long- wave fluctuations are reduced). This corresponds to Debye shielding in plasma physics in 
which each particle is surrounded by a cloud of opposite sign particles. When /3 < 0, the spec- 
trum is enhanced at small k (i.e. long- wave fluctuations are increased). This corresponds to a 
form of anti-shielding in which each particle is surrounded by a cloud of similar particles. This 
is similar to the Jeans clustering in astrophysics. The second term in the correlation function 
( 18 ip arises from the vortex cloud which, on average, surrounds each particle. The radius of 
the cloud is L ~ k^ l {l — /3/A)) -1 ^ 2 and the effective number of vortices contained in the cloud 
(obtained by integrating (151]) ) is N e ff ~ (1 — /? / For (3 > 0, the cloud is of opposite 
circulation to the vortex and shields its effect. For (3 < 0, the cloud is of same circulation as the 
vortex and enhances its effect. The effective number of vortices in the cloud increases as the 
inverse temperature decreases and it diverges for f3 = (3q. This corresponds to the formation of 
large "clumps" or "clusters" of vortices. For (3 < Pq, the statistically homogeneous distribution 
of vortices becomes unstable and forms clusters. When the system becomes spatially inhomo- 
geneous, the dynamics of clustering can be studied with the mean field equation (1581) coupled 
to the screened Poisson equation Aip — k^tp = —uj. 

2.8 The different regimes and the orbit-averaged-Fokker-Planck equa- 
tion 

Using the previous kinetic theory, we can qualitatively discuss the different regimes of the 
dynamics of 2D Brownian vortices depending on the values of N, v and \i = (3v (a similar 
discussion applies to self-gravitating Brownian particles jlQ] and to the BMF model [58] ) . For 
v = ji = and N — > +oo, Eq. (|2"T|) reduces to the 2D Euler equation. Starting from an unstable 
initial condition, the 2D Euler equation can undergo a process of violent relaxation and form 
a quasi- stationary state (QSS) on a short dynamical timescale ~ to [591 ED]- Then, finite N 
effects are expected to drive the "collisional" relaxation of the point vortex gas towards the 
microcanonical Boltzmann distribution [with temperature /3(E)] on a longer timescale t re i ax (N) 
whose scaling with N is not yet known (see discussion in [2^1 ED])- For v > 0, \i ^ and N 3> 1, 
the mean field Fokker-Planck equation (127]) drives the system towards the canonical Boltzmann 
distribution (1281) [with bath temperature j3\ on a diffusive timescale tdiff ~ R 2 jv independent 
on N. We can distinguish three cases: (i) When to <C t re i ax (N) <C tdiff (corresponding to 
N finite and u, /j, —* 0), the system undergoes violent relaxation towards an Euler QSS, then 
evolves towards a microcanonical quasi stationary state [58] due to "collisions" (i.e. correla- 
tions due to finite N effects) and finally relaxes towards the canonical equilibrium state, (ii) 
When t D t di ff t re i ax (N) (corresponding to N — ► oo and relatively small viscosity), the 
system undergoes violent relaxation towards an Euler QSS, then relaxes towards the canon- 
ical equilibrium state without forming a microcanonical quasi stationary state, (iii) When 
tdiff ^ tn <C t re i ax (N) (corresponding to large viscosity v,(/,—> +oo), the system undergoes a 
diffusive relaxation towards the canonical equilibrium state without forming a QSS. 

Let us consider more specifically the second case. We assume that N — » +oo so that 
"collisional" effects are completely negligible. We also assume that v — ► and \i — > (while 
(3 ~ 1). To a first approximation, the kinetic equation (l27jl reduces to the 2D Euler equation. 
We assume that the system has reached a QSS as a result of violent relaxation. This is a stable 
stationary solution of the 2D Euler equation of the form u = f{ip), which depends only on 
the stream function. It is reached on a few dynamical times tp. If v is finite, the system will 
evolve on a longer timescale due to the terms of drift and diffusion present in the Fokker-Planck 
current. If v and \i are sufficiently small, the solution uj(r,t) will remain close, at any time, to 
a steady solution of the 2D Euler equation of the form uj = oo(ip,t) that slowly evolves in time 
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due to the effects of diffusion and drift. Noting that 



we can rewrite the kinetic equation (1271) in the form 

du dip du 
dt dt dip 



d , , du dip du 

-u^M),t) = - + -- ; 



Q(u) 



(83) 



(84) 



where Q(u) = — V ■ J is the Fokker-Planck term. Using the fact that u only depends on the 
stream function in a first approximation, we have 



Q(u) = uV ■ (Vu + (3^uVip) = vV 



du 
dip 



Pju \7ip 



v + (3ju ) Aip + uVip ■ V 



uAip 



\dip 
du 
dip 



du 
dip 



P^u 



+ ) + u(vipy 



d f du 



dip \dip 



+ P'ju 



(85) 



where Aip = k\ip — u for the QG model (Aip = —u for the ordinary point vortex model). The 
next step is to get rid of the dependence on r in Eq. ([85"]) by averaging over the orbits. Then, 
the dynamical equation for u(ip,t) becomes 



du I dip\ du 
~dt + \~dt/ Ihp 



( du 



\dip 



d f du 



u(k 2 ip -u)[— + f3ju)+v ((ViP) 2 ) —l— + p 7U ), (86) 



dip \dip 



where (X) denotes the average of X(r, t) on the elementary surface of the plane between ip and 
ip + dip. Defining 



x y>,t) = <(V^) 2 > 



Wr,t)<W V ^) 2rfr _ JS(iP-iP(r,t))(ViP) 2 dr 



il><ip(r,t)<il>+dil> 



dr 



f6(iP-iP(r,t))di 



»7) 



and using (dip/dt) = d{ip)/dt = 0, we finally obtain the orbit-averaged- Fokker-Planck equation 



d_ 

dt 



( du 



u(ip, t) = u(k 2 ip - u) — + + u X (ip, t)— — + pyw , 



\dip 



d f du 



dip \dip 



Aip = k 2 ip - u(ip,t). 
The steady solution of this equation is the Boltzmann distribution (1281) . 



(89) 



3 The two-species system 
3.1 The stochastic equations 

For a multi-species system of point vortices, the stochastic equations consistent with the Gibbs 
canonical distribution at statistical equilibrium are 

— - = z x ViH- uViH + v^Ri(t), (90) 

dt ji 
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where H = J2i<j1i1j u ( v i^ r j) + (1/2) J2iLi lf v { r i-> r «) is the Hamiltonian. The Fokker-Planck 
equation determining the evolution of the iV-body distribution function Pjv(ri, rjv, t) asso- 
ciated with the stochastic equations ([90]) reads 

dPN ^STv 9Pn ST 8 ( 9Pn v dH \ rem 

^r + E v -^- = E5r--^^- + ^^-J' ( 91 ) 

i=l 1 i=l 1 ^ 1 1 / 

where Vj = — (l/7j)z x OH/dvi is the total advective velocity of point vortex i. For z/ 7^ and 
fi 7^ 0, the stationary solution of the Fokker-Planck equation fl9lT) is the canonical distribution 



, 1 

N 



Z e-? H , (92) 



provided that the mobility and the viscosity are related to each other through the Einstein 
relation ([8]). The free energy of the iV-body system is given by Eq. (Q and we have the same 
general properties as for the single species system (see Sec. 12. 2p . If we restrict ourselves to 
positive temperatures (i.e fi > 0), the stochastic equations (|9"U|) without the rotational term 
provide a microscopic description of the Debye-Huckel model of electrolytes [31]. At positive 
temperatures, the full equations (1901 also provide a microscopic description of a dissipative 
guiding center plasma under a strong magnetic field, where like-sign charges repel each other 
(see Appendix |AJ . Here, we shall be more general and consider both positive and negative 
temperatures, i.e repulsive or attractive interactions between like-sign particles. 

To be specific, let us consider a globally neutral two-species system made of N/2 point 
vortices of circulation +7 and N/2 point vortices of circulation —7. For box-confined vortices 
interacting through the Newtonian (or Coulombian) potential fl56|) . it can be shown that sta- 
tistical equilibrium states exist ^ in the canonical ensemble (i.e. the partition function is finite) 
iff (see [2S1 EE] and references therein): 



ft" - < /» < K = 03) 



At positive temperatures (3 > (implying ji > 0), the interaction is "repulsive": the point 
vortices of the same sign have the tendency to repel each other while the point vortices of 
opposite sign have the tendency to attract each other. At large temperatures (small (3) a 
vortex of a given circulation is surrounded by a cloud of vortices with opposite circulation 
which screen the interaction. This is like the Debye shielding in plasma physics. At small 
temperatures (large j3) the vortices have the tendency to form pairs or microscopic dipoles 
(+, — ) similar to "atoms" (+e, — e) in plasma physics. In this regime, the flow consists in a 
spatially homogeneous distribution of N/2 individual dipoles. At the critical temperature 
the most probable state is a gas of N/2 singular microscopic dipoles (+, — ) where the vortices 
of each pair have fallen on each other. This leads to N/2 Dirac peaks made of two vortices 
of opposite circulation. Thus, the value of the critical temperature (3~£ can be understood by 
considering the statistical mechanics of only two vortices (one with positive circulation +7 and 
one with negative circulation —7) and determining at which temperature the partition function 
ceases to be normalizable. The critical temperature /3£ can thus be obtained from Eq. (1571) 



8 The equation of state is known exactly for a 2D point vortex gas (see, e.g., [35]). Note that the pressure 
must be positive at positive temperatures (3 > and negative at negative temperatures (3 < 0. This determines a 
range of forbidden temperatures where the pressure has not the right sign: (3 > (3^ = 8ir/j 2 for neutral systems 
and (3 < (3~ — — &n/[(N — 1)7 2 ] for a single species system of N point vortices. According to inequalities (|57|) 
and (|93p . the partition function diverges strictly before entering this range of temperatures, so these states are 
not accessible anyway (see, e.g., [38]). 
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by taking N = 2 and changing the sign (since the "attraction" of like-sign vortices at (3 < 
corresponds to the "attraction" of opposite-sign vortices at (3 > 0). For (3 > there is 
no statistical equilibrium state anymore and this regime can be studied dynamically with the 
stochastic model (!90l) . At negative temperatures f3 < (implying \x < 0), the interaction is 
"attractive" : the point vortices of the same sign have the tendency to attract each other while 
the point vortices of opposite sign have the tendency to repel each other. At sufficiently small 
temperatures (sufficiently negative j3), the vortices have the tendency to form two macroscopic 
clusters: one cluster of N/2 positive vortices and one cluster of N/2 negative vortices. In this 
regime, the flow consists in a spatially inhomogeneous distribution of point vortices forming 
a macroscopic dipole. At the critical temperature f3~ , the most probable state is a singular 
macroscopic dipole (y+> y~ ) wriere the vortices of the same species have fallen on each other, 
creating two Dirac peaks made of N/2 vortices. The value of the critical temperature j3~ 
can be understood by considering the statistical mechanics of a single species system of N/2 
vortices and determining at which temperature the partition function ceases to be normalizable. 
The critical temperature f3~ can thus be obtained from Eq. (157)) by making the substitution 
iV — > N/2. For (3 < (3~ there is no statistical equilibrium state anymore and this regime can 
be studied dynamically with the stochastic model (19"U)) . 



3.2 The BBGKY-like hierarchy 

For sake of generality, we first consider a two-species system made of N + = Nn + point vortices 
of circulation 7+ (labeled from i — 1 to n + N) and N_ = Nn_ point vortices of circulation 
7_ (labeled from i = n + N + 1 to N). We use the same notations and presentation as in the 
equilibrium study of Pointin & Lundgren [T5]. The one-body distributions are defined by 

P+(r u t)= [ P N (r 1 ,...,r N ,t)dr 2 ...dr N , (94) 



P 1 (r N ,t) = j P N (r 1 ,...,r N ,t)dr 1 ...dr N _ 1 . (95) 

The two-body distributions P 2 ++ (r l7 r 2 , t), P 2 +_ (r 1 , r 2 , t) = P 2 ~ + (r 2 , r 1; t), P 2 (ri,r 2 ,£), and 
the higher-body distributions are defined similarly. From the Fokker-Planck equation ()9ip . it 
is easy to derive the complete BBGKY-like hierarchy of the two-species point vortex gas. The 
general equations are 







i=l 



d ( dP 
1 v 



dvi 



HP} 



p T ,n' \ ^ 2 



dri 



p , n 



E 

k=p+l 



7+7- 



9'Uki ( p\ 2 f p (p+l) + ,n- 9'u kl 



n 



+»N{n---) 1+1 - \P^[ 



,p+,(n+l)- 



d'u ki 
<9r ?: 



& 



E 



-)> (96) 



i=p+l 



where j = p + n. The first equation of the BBGKY-like hierarchy for P x + is 



dPf 
~dt 



v 



_d_ 
dri 



+/3A^n_7 + 7_ 



dvi 



& 



P 2 +-(l,2Mr 2 + -/3P+7^^(l,i; 



(97) 
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The equation for P 1 is obtained by interchanging + and — . The second equations of the 
BBGKY-like hierarchy for P 2 ++ an d P£~ are 



dt 



v 



d_ 
dri 



/3 7 ?P 2 ++(1,2)^ + /37V 



+/3iV?7,_7 + 7_ 



^P 3 ++ -(l, 2, 3)dr 3 + ^^^(1, 1) 



(—1)^/ ^-(1,2,3)^ 

+ (1-2), (9* 



<9ri 



+ ^-Pt-Q, 2)^f + BN („ + - I) 7 * / ^fP 3 +-+0, 2, 3 )<ir3 



W [n- - i) 7+7- y 



^P 3 + --(l,2,3Mr 3 + ^P + -7 + ^(l,l) 



+ (!<-> 2). (99) 



The equations for P 2 and P 2 + are obtained by interchanging + and — . Next, inserting the 
Mayer decompositions 



P 2 ++ (l,2)=P 1 + (l)P 1 + (2)+P 2 ++ (l,2), 



(100) 



P 2 + -(l,2)=P 1 + (l)Pr(2)+P 2 + -(l,2), 
in Eq. (|97|) . we obtain the exact equation 



dP? 



_d_ 
dri 



+(3N\n + --) 



&U12 



dri 

+/3iV?7,_7 + 7 



P 2 ++ (l, 2)rfr 2 + (3Nn. 1+1 .P+{\) 



&U12 

dri 



Pf (2)rfr 2 



^P 2 -(l,2)rfr 2 + ^P + 7 + ^(M) 



Inserting the Mayer decompositions 



P 3 + ++ (l, 2, 3) = P+(l)P+(2)P+(3) + P+(l)P 2 ++ (2, 3) + P 1 + (2)P 2 ++ (1, 3) 

+P 1 + (3)P 2 ++ (1,2) + P 3 ' +++ (1,2,3), 



(101) 



(102) 



(103) 



P 3 + + -(l, 2, 3) = P 1 + (l)P 1 + (2)Pf(3) + P 1 + (l)P 2 + -(2, 3) + Pf{2)P'+-{l, 3) 

+Pr(3)P 2 ++ (l,2)+P 3 ' ++ -(l,2,3), 



(104) 



P 3 +-+(l, 2, 3) = P+(l)Pf(2)P 1 + (3) + Pi-(l)P 2 - + (2, 3) + Pf (2)P^ + (1, 3) 

+P 1 + (3)P 2 + -(1,2)+P 3 + - + (1,2,3), 



(105) 



P 3 +--(l, 2, 3) = P+(l)Pf (2)Pf(3) + P 1 + (l)P 2 "-(2, 3) + P 1 -(2)P 2 +-(1, 3) 

+Pr(3)P 2 + -(l,2)+P 3 + -"(l,2,3), 



(106) 
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in Eqs. 



, and using Eq. (j!02p to simplify some terms, we obtain the exact equations 
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(107) 
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•(l,2)/^P 1 -(3)dr 8 
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+3 7+7 _Pf (2) | ^f^2 + -(l, 3)rfr 3 + (n_ - 1) 7+7 _P 2 - 

r* 1 
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To order 1/N in the proper thermodynamic limit defined in Sec. 12.21 the foregoing equations 
reduce to 



+ PN [n 



1)^(1)/^)*, 
+/?iVn +7 * | ^Hp; ++ (l,2)rfr 2 + /3iVn_ 7+7 _P 1 + (l) | ^Pf(2)dr 2 
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dP, 



dt 



_d_ 

dri 



+ /37 + 7-P 1 + (l)Pf(2)^ + /5iVn +7 ^(l) f ^Hp;-+(2,3)dr 3 



Or i 



Or i 



Or i 



+PNn + T 2 + P f +-(l,2) |^P 1 + (3)dr 3 -/3 7+ 7-A + (l)^r(2) J ^fPf(3)dr 3 



+/3iVn_7 + 7^i + (l) y^^-(2,3)dr 3 + /3JV n _7+7-^ + -(l,2) | ^fPf(3)eZr 3 



+ (!<-> 2). 



(Ill) 



Before considering simplified forms of the above equations, we shall derive useful expressions 
of the average energy. 

3.3 The average energy 

The average energy 

f 1 N f 

E = (H) = / liljU ij P N dr 1 ...dr N + - ^ / l 2 i v ii P N dr l ...dv N , (112) 

can be written 

E = ]p 2 + N + (N + -l) j P 2 ++ (l,2)n 12 rf ri rfr 2 + i 7 !iV_(iV_-l) ^ P" (1, 2)u l2 dv 1 dv 2 

+7V + iV_7 + 7_y P 2 + "(l,2) Ml2 rfr lC /r 2 + l | [^ +7 JP+(1) + JV_7iPf(l)] l)dn. (113) 
Substituting the Mayer decomposition in the previous equation, we obtain 

P+(l)P+(2) +P2++(1,2)] uudndr, 
Pf (l)Pf (2) + P 2 — (1,2)] u 13 d ri dr 2 
P+(l)Pf(2) + P 2 ' + ~(l,2)] u^dndra 

y K7+A + (l) + ™~7 2 iT(l)] u(l, (114) 

Introducing the average vorticity 

u(r, t) = Nn+7+P+(r, t) + JVn_7_Pf(r, i), (115) 
and the corresponding stream function ( 123|) . the average energy can be put in the form 

n+7+P 2 ++ (l, 2) + n_7!P^--(l, 2)] u^dv^ 

7>+P 2 ++ (l, 2) + 2n + n_ 1+1 _P' 2 + ~(l, 2) + 7^^— (1, 2)] u 12 dr 1 dr 2 

J P+{l)P+{2)u l2 dT l dr 2 -]plNn. J P{{l)P{{2)u 12 dv 1 dv 2 
N f 

+- J [n +7 ^P 1 + (l) + n_7^Pr(l)] u(l, l)dn. (116) 
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This equation is exact for any N. At the order 1/N in the proper thermodynamic limit defined 
in Sec. 12.21 it reduces to 



If N f 

E = - ui>dr+- / [n +1 2 + P+(l] 



- n 



_ 7 !Pr(i)] u(i,i)dn 



4^ 



7+<P2 ++ (l, 2) + 2n + n„ 7+7 „P2 + "(l, 2) + TinlPg— (1, 2)1 u 12 dndr 2 



— 7+JVra. 



. y P+(l)P+(2)« 12 dr 1 dr 2 - ^7-^- / P{{l)P{{2)u 12 dr 1 dr 2 . (117) 



3.4 The mean field approximation 

In the limit N — > +oo, the first equation fllQ9p of the BBGKY-like hierarchy reads 
dPt 



Of 



"^{^f + /?7+Pl+(1) / [Nn + l + Ptm+Nn^P { {2)] dr 2 }. 
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Introducing the average vorticity (11151) and the corresponding stream function (123]) . we can 
rewrite Eq. (1 11 8 [) in the form 



<9P+ 



+ u ■ VP X + = z^V • (VP^ + /37+P+V^) 



119) 



where u = — z x is the mean advective velocity. The equation for Pf is obtained by 
interchanging + and — . Then, defining u + = Nn + '~f + Pi and u;_ = Nn-j-P^ so that = 
cj + + cj_, and considering a potential of interaction of the form fl56l) . we obtain the coupled 
system of equations 



duj + 



+ u • Vw+ = i/V ■ (Vcu+ + /37+cj+V^) = -V • J. 



(120) 



<9u;_ 



+ u • Vc^_ = uV ■ (Vcu_ + /37_cj„V^) = -V • J. 



I2r 



Aip = u + + U- 



122) 



with the normalization J u + dr = Nn + '~f + = T + and J oj_dr = iVn_7_ = T_. These mean field 
Fokker-Planck equations conserve the circulations (or vortex numbers) of each species. The 
steady solutions are the Boltzmann distributions 



A e -/37+V> 



A_e~^. 



(123) 



They can also be obtained from the equilibrium BBGKY hierarchy in the limit iV — > +oo. Sub- 
stituting Eq. (I123p in Eq. HI 22)) . we obtain a mean field equation determining the equilibrium 
stream function. For a globally neutral system with 7 + = — 7_ = 7 and n + = n_ = 1/2, this 
equation reduces, in the symmetric case A + = —A_, to the celebrated sinh-Poisson equation 
Aip = 2 A smh((3 , yip) studied in [151 [61] (note, however, that there is no fundamental reason why 
A + = —A-] this corresponds to very particular initial conditions). At positive temperatures, 
and for a globally neutral system, the kinetic equations fll20l) - fll22l) . without the advective 
term, are isomorphic to the Debye-Hiickel model of electrolytes where like-sign charges repel 
each other (see Eq. (3) of [34]). Therefore, for (3 > 0, our study provides a kinetic derivation, 
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from the microscopic stochastic process f[§Ul) . of the mean field kinetic equations introduced by 
Debye & Hiickel [M]. At positive temperatures, these equations lead at statistical equilibrium 
to a spatially uniform state with zero net charge. The novelty of our model is to allow also for 
the consideration of negative temperatures where like-sign vortices attract each other. For suf- 
ficiently negative inverse temperatures, these equations lead to a macroscopic order, typically 
the formation of a large-scale dipole |62j . 

The Boltzmann entropy of the two-species point vortex gas 



S = -N J [n+Pf(l) hiP 1 +(l) + n_Pf (l)lnPf(l)] dx u 



(124) 



can be rewritten (up to an additive constant): 



S = - / — In -i + — In — dr. (125) 
J V7+ 7+ 7- 7-/ 

This expression of the Boltzmann entropy can also be obtained from a classical combinatorial 
analysis [30J. On the other hand, for iV — > +oo, the mean field energy is 

E = - I u^dv. (126) 



2 . 

Then, the mean field free energy (more properly the Massieu function) reads 

J = S-/3E, (127) 

where S and E are given by Eqs. (11251) and (11261) . The mean field free energy (11271) can be 
obtained from Eq. (Q in the limit N — > +oo. It is straightforward to establish that 



J= I I — + — \dr > 0. (128) 

This is the appropriate form of iJ-theorem in the canonical ensemble: J > and J = iff uj± 
are the mean field Boltzmann distributions (11231) . The free energy J[a> + ,u;_] is the Lyapunov 
functional of the Fokker-Planck equations (I120l) - (ll22p . The Boltzmann distributions (I123p are 
critical points of J at fixed r + and T_ and they are linearly dynamically stable iff they are (local) 
maxima of J. If J is bounded from above, we conclude from Lyapunov's direct method, that 
the mean field Fokker-Planck equations (1 1 20[) - (1 1221) will reach, for t — > +oo, a (local) maximum 
of J at fixed circulations T + and T_. If several local maxima exist (like dipoles and bars in a 
periodic domain [63]), the selection of the maximum will depend on a complicated notion of 
basin of attraction. If there is no maximum of free energy at fixed circulations, the system will 
have a peculiar behavior and will generate singularities. This is the case in particular for (3 < (3~ 
in the situation described in Sec. 13.11 In that case, the partition function diverges and the mean 
field free energy has no global maximum. This situation can be studied dynamically by solving 
the mean field Fokker-Planck equations (1 1 20 p - (1 1 2 21) . Preliminary numerical simulations of 
Eqs. (1 1 20 j) - (1 1 221) show the formation of a singular macroscopic dipole^. The results of these 
simulations will be reported elsewhere [62]. We note that the relaxation equations (I120p - (I122I) 
can be obtained from a maximum entropy production principle (MEPP) by maximizing the 
production of free energy J at fixed circulations T± (and additional physical constraints) [51j. 



9 Note that a similar growing condensate (dipole) has been studied recently in the context of forced 2D 
turbulence [64] . Our model of Brownian vortices is physically different from [64] but the connexion between the 
two studies may be noted. 
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It is also straightforward to generalize the present results to a multi-species system of Brow- 
nian vortices. The general mean field Fokker-Planck equations are given by 



OuJa 

dt 



(129) 



where ip is given by Eq. (1231) with uj = ^2 a uj a where u a (r,t) = N a '-f a P^ a \r,t) is the ensemble 
averaged vorticity (proportional to the density) of species a. The free energy (Massieu function) 
is J = S - 13 E where S = - ]T a / ^ hi ^ dr and E is given by Eq. (TO . The ^-theorem 
reads 

j = v V / (Voo a + (3^ a uj a Vtp) 2 dr. (130) 

The relaxation equations (11291) have the following properties: (i) T a = J u a dr are conserved 
(ii) J ~> (iii) J = iff u a = A a e~^ la ^ (iv) a steady state is linearly stable iff it is a (local) 
maximum of J at fixed T a . Thus, if J is bounded from above, these equations will relax towards 
the canonical statistical equilibrium state of the multi-species point vortex gas. We can also 
extend the approach of Sec. 12.51 to obtain the equations satisfied by the exact vorticity field of 
each species. They are given by 



~dt 



u • Vu a = z/V • (Vu a + (]^ a UJ a Vlp) + V ■ (y/2jsy c 



UJ, 



'131] 



where ip is given by Eq. PSj) with 10 = ^ a w a where u a (r,t) = 7«X]ieY Q — r «(^)) i s the 
exact vorticity field of species a (here, X a denotes the ensemble of point vortices of species 
a). Using the same arguments as those given at the end of Sec. 12. 5[ we can introduce a 
coarse-grained vorticity field that smoothes out the Dirac distributions while keeping track of 
stochastic effects. The evolution of the coarse-grained vorticity is then given by an equation of 
the form (I13ip where ip is given by Eq. (1231) with u = Yla^a where u Q (r,t) denotes now the 
coarse-grained vorticity of species a. These stochastic kinetic equations, involving a noise term, 
could be used to study dynamical phase transitions between different maxima of free energy, 
like "bars" and "dipoles" in periodic domains. This will be considered in future contributions. 



3.5 The two-body correlation function 

We now consider a globally neutral system such that 7+ = +7, 7_ = —7 and n + = n_ = 1/2. 
We also assume that the distribution of point vortices is spatially homogeneous (to leading 
order) so that P 1 ± (r, t) = P + P ± (r,t) where is of order 1/N. In that case, the second 
equations (jliop - fllllj) of the BBGKY hierarchy at the order 1/N reduce to the forms 



dP' 2 ++ 

dt 

+ l -f3N 1 2 P a 



d (dP 



V 



dri 
d'u i3 
dr 1 



dr 1 



P; ++ (2,3)-P 2 + -(2,3) 



dr 1 



dr 3 



d'u 



13 



dr- ? 



dr 1 

(1^2) 



(132) 



dP' 



d (dP 



dt 

+^iV 7 2 Po 



v 



dv\ 
d'ui 3 
dri 



dri 



OTi 



PT+(2,3)-P;-(2,3) 



dr 3 



d'u 13 



dr ? 



9ri 

(1-2), 



(133) 
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and the energy of interaction is 



iv 2 7 2 



E 



P;++(1 ) 2)-2P 2 +-(1,2) + P 2 --(1,2) 



NY 



J Pfaudv^ + ^N^Po J u(l,l)dn. 



(134) 



Considering solutions such that (note the presence of the term 1/N in the definition of h): 



P' 2 ++ {1,2) = P;-(l,2) = -P; + -(l,2) = -P^ + (l,2) = P 2 



h(r 1 -r 2 ,t) + — 



(135) 



we find that the two-body correlation function h(r\ — r2,t) satisfies an equation of the form 

This returns Eq. (|35|) for a one component system with a neutralizing background. The 
correlational energy is 



E = E- ^N^Pq 



2„,2 



P 2 / h(l,2)ui 2 dr 1 dr 2 , 



(137) 



which coincides with Eq. fl65|) . Finally, the first equation of the BBGKY-like hierarchy, written 
at the order 1/N, gives the evolution of the function P + (r,t) in the form 



dP + d 
dt drx 



8P+ 
dr 1 



-BN^Po 



dr 1 



P + (2)-p-(2) + 2/i(l,2)P 



dr? 



-M 2 Pl 1 d ' Ul2 "- 



(138) 



and an equivalent equation for P~{v,t). Since we recover the same results as in Sec. 
this two-species model justifies, on a firmer basis, the results obtained in Sees. 12.61 and 12.71 
based on a "Jeans-like swindle". We note that the correlation function fl6"3"l) and the caloric 
curve (166]) are defined for all positive temperatures. In particular, these formulae seem to be 
valid above the critical inverse temperature 8^. This is of course incorrect and this clearly 
shows that the Debye-Hiickel theory is not valid for low temperatures (note that the random 
phase approximation in [T7] also breaks down at low temperatures). In particular, it cannot 
account for the formation of tightly bound dipoles or pairs (+, — ) as 3 — > 0^. To arrive 
at the expressions (1631) and (1661) . we have assumed that the three-body correlation function 
P3 could be neglected because it is of order 1/N 2 . This is correct "on average" but this 
is not correct at any scale, in particular for small separations. Therefore, it is necessary to 
take into account higher order correlation functions in the kinetic theory (and go beyond the 
Debye-Hiickel approximation) if we want to reproduce the pair formation at low temperatures. 
In this respect, we note the important paper of Samaj & Travencc [65J who obtained the 
exact thermodynamic parameters of the 2D two- components plasma for any value of accessible 
positive temperatures < Q < This study clearly shows that the Debye-Hiickel theory is 
only valid at large temperatures. 
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4 Conclusion 



In this paper (see also [32J ) we have introduced a new model of random walkers in long-range 
interactions that we called "Brownian vortices" . By definition, their dynamics is governed by 
stochastic equations of the form ( 1901) . This dynamical model is associated with the canonical 
ensemble and it leads to the Gibbs canonical distribution at statistical equilibrium. In the case 
of Brownian vortices, we can develop a kinetic theory and analyze all stages of the dynamics. 
In the mean-field limit, we find that the smooth vorticity field is governed by a parabolic 
drift-diffusion equation (J27j) coupled to an equation of the form (1231 . which is often an elliptic 
equation like the Poisson equation. This type of equations is actively studied at present by 
physicists and applied mathematicians in relation to self-gravitating Brownian particles and 
biological populations experiencing chemotaxis [66] • We think that the present model naturally 
enters in this line of investigations. In particular, an interesting novelty of this model is to 
allow for particles (vortices) with positive and negative circulations (while the mass of material 
particles -like stars or bacteria- is always positive) leading to a wider richness of phenomena. 

Another interest of the point vortex model is to allow for both positive and negative tem- 
peratures (see Table [2]). At negative temperatures, the interaction between like-sign vortices is 
"attractive" (like gravity) and can create a macroscopic order characterized by the nonunifor- 
mity of the spatial distribution of the vortices (typically a monopole or a dipole). The statistical 
mechanics of point vortices at negative temperatures shares a lot of analogies with the statistical 
mechanics of stellar systems in astrophysics [53]. It is, furthermore, richer because we can have 
positive and negative circulations resulting in the separation of a cluster of positive circulations 
and a cluster of negative circulation (dipole) while gravitational systems typically organize in 
a single cluster. We have also mentioned the existence of a critical inverse temperature (3~ 
below which the partition function diverges and the free energy is not bounded anymore. In 
that case, we expect the system of point vortices to collapse. It forms a singular monopole A7 
(Dirac peak) in the single-species case and a singular dipole (y7, — y7) in the two-species case. 
At positive temperatures, the interaction between like-sign vortices is "repulsive" (like in elec- 
trostatics). For a neutral system, this leads to a spatially homogeneous state with correlations 
which behaves therefore like a 2D Coulombian plasma. We have mentioned the existence of a 
critical inverse temperature fl£ above which the partition function diverges and the free energy 
is not bounded anymore. In that case, we expect the system of point vortices to collapse. It 
forms a gas of N/2 singular dipoles (+7, —7) similar to atoms (+e, — e) in a 2D plasma. 

For inhomogeneous systems, the mean field energy (j!26p scales like N 2 , i.e. the normalized 
energy e = E/{N 2 ^ 2 ) is of order unity, corresponding to a nonextensive scaling. This regime 
of negative temperatures, leading to a macroscopic order (like for self-gravitating systems), 
has been studied by Onsager [24] and others PSI HSl D21 UHl HS1 1213 I2H 1221 12B] at statistical 
equilibrium. On the other hand, at positive temperatures, the interaction between like-sign 
vortices is "repulsive" (like for electric charges) and the physics is very similar to that of a 
two-dimensional plasma. In that case, we recover the model of Debye-Huckel electrolytes [54"] . 
In the neutral case, the distribution of vortices is spatially homogeneous due to Debye shielding. 
For homogeneous systems, the correlational energy (11371) scales like N (since P' 2 ~ 1/A), so it 
has an extensive scaling. This regime of positive temperatures, where the usual thermodynamic 
limit applies, has been studied by Ruelle [67] and others [TTJ [68] at statistical equilibrium. It was 
believed for some time that these two approaches (Onsager and Ruelle) were in contradiction. 
In fact, they correspond to different definitions of the thermodynamic limit describing different 
situations. The present model of Brownian vortices provides an out-of-equilibrium model that 
is fully consistent with the results of statistical equilibrium in the canonical ensemble at both 
positive and negative temperatures in the two situations described above. 
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The perspectives of this work are the following. On a theoretical point of view, it would be 
important to study numerically, analytically and with the methods of applied mathematics the 
mean field Fokker-Planck equations fl 1 2 [) - (11 2 2 j) . This has largely been done already in the single 
species case since the corresponding equations (|58l) - (l59l) are isomorphic to the Smoluchowski- 
Poisson system and Keller-Segel model [32] • However, the two species mean field Fokker-Planck 
equations (I120p - (I122I) have not been studied extensively so far (especially in the attractive case 
valid in the regime of negative temperatures). In particular, the formation of the singular 
macroscopic dipole for (3 < (3~ is of particular theoretical interest [62J. The study of the N- 
body stochastic equations (j2J) or (1901 and the study of the stochastic kinetic equation (1551) or 
( 11311) are also of interest. When the free energy J has several local maxima at fixed circulation 
(metastable states), representing for example monopoles, dipoles, bars or jets like in [63], the 
stochastic term (noise) in Eq. (|55|) or (11311) can induce dynamical phase transitions from one 
state to the other. It would be interesting to study this phenomenon in detail. On the other 
hand, on a physical point of view, it would be interesting to know if we can experimentally 
realize a system of Brownian vortices described by Eqs. (I90p . It would be interesting to devise 
experiments (not necessarily in the context of fluid mechanics) where the effective coupling of 
the particles with a thermal bath can be modeled by stochastic processes of the form (1901) . In 
the repulsive case, these equations can describe a dissipative plasma (see [31] and Appendix 
|A|) . In the attractive case, these equations could describe biological systems where two species 
of particles organize in two distinct clusters, i.e the particles of the same species attract each 
other and the particles of different species repel each other. 

Let us finally comment on the inequivalence of microcanonical and canonical ensembles for 
the point vortex system (see also Appendix [C]). For ordinary point vortices whose dynamics is 
governed by Hamiltonian equations, the statistical equilibrium state is obtained by maximiz- 
ing the Boltzmann entropy Sb at fixed energy E and circulation V (microcanonical ensemble). 
For Brownian point vortices whose dynamics is governed by Langevin equations, the statistical 
equilibrium state is obtained by maximizing the Boltzmann free energy J# = Sb — /3E at fixed 
circulation V (canonical ensemble). Canonical stability implies microcanonical stability but the 
reciprocal is wrong. This means that some states (E,/3) that are accessible in the microcanon- 
ical ensemble (maxima of entropy S) may be inaccessible in the canonical ensemble (they are 
not maxima of free energy J). When this happens, we speak of ensemble inequivalence. At 
equilibrium, ensembles are equivalent for vortices in a disk (or other simple domains) when the 
Hamiltonian is postulated to be the only constraint [221 EH]- However, they can be inequiv- 
alent in other kinds of domains [69] or also in a disk when the conservation of the angular 
momentum is taken into account [TQl [23] . Out-of-equilibrium, the kinetic equations describing 
Hamiltonian (microcanonical) and Brownian (canonical) point vortices are very different in 
any case. For Hamiltonian vortices, when N — > +oo, the kinetic equations reduce to the 2D 
Euler equation and finite N effects are necessary to get the convergence to the microcanonical 
Boltzmann distribution with temperature (3(E) [26J. The relaxation time is larger than Nto 
and its precise scaling with N is not known (we are not even sure that the system will reach 
thermal equilibrium). For Brownian vortices, when N — > +oo, the kinetic equations reduce 
to the mean-field Fokker-Planck equations (|27p or (I120l) - (ll22p . They converge towards the 
canonical Boltzmann distribution with the bath temperature (3 on a diffusive timescale L 2 /v 
independent on N. Therefore, although the equilibrium states correspond to mean field Boltz- 
mann distributions in the two ensembles, the convergence to these steady states is radically 
different in the Hamiltonian and Brownian descriptions. Thus, our study explicitly shows the 
fundamental analogies and differences that exist between a microcanonical and a canonical de- 
scription of point vortices. These analogies and differences are summarized in Table [T] where 
we compare the kinetic equations of Hamiltonian and Brownian vortices (a completely similar 
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comparison could be made between Hamiltonian and Brownian systems of material particles 
with long-range interactions [U El [6j [7J [8] ) . 



A Dissipative guiding center plasma 

The equations of motion of a system of charges in a vertical magnetic field B = Bz are 

dvi V; x B 

m— = e i E i + e i — - — , (139) 

where Ej = — Vj$ = ^- X)jyi(VO e j lnx^- is the electric field created self-consistently by the 
charges [TS]. For large fields B ^> 1, we can neglect the inertial term dvi/dt ~ and make the 
guiding center approximation 

v, ~ ~B x Ei. (140) 



The equations of motion are then similar to those of a Hamiltonian 2D point vortex gas [15] . For 
a conservative system (microcanonical ensemble) both positive and negative temperature states 
are possible [21]. Let us now consider a dissipative magnetized plasma so that the equations of 
motion can be modeled by the stochastic process 

£=£(* + ^)-«*. + V5«R,( t ), (141) 

where £ > is an effective friction coefficient and Rj(i) a random force with (Rj(t)) = 
and (R?(t)Rj(t')) = 8ij8 a p8{t — t'). The Einstein relation reads £ = Dm/3, implying that 
the temperature is necessarily positive for this model. If we consider both a large field limit 
B — > +oo and a strong friction limit £ — > +oo in such a way that £ ~ eB/(mc), we can neglect 
the inertial term in Eq. (11411) . Then, after simple manipulations, we obtain the stochastic 
equation 

m2C )e) v, = ^E, - ^B x E + ^ V2DR, - ^V2DB x R, (142) 



efB^ J 1 aB 2 1 B 2 e 2 B 2 1 ei B* 

Defining Qj = a«Rj — B x Rj with a« = £mc/e.j, we have (Qi(£)) = and (Qf(t)Q?(t')) 



2 



{ai + I)8ij8 a p8(t — t'). Therefore, the stochastic equation (11421) can be written in the form of 
Eq. (1901) provided that the parameters are properly reinterpreted. If we consider a dissipative 
plasma, without imposed magnetic field (B = 0), the stochastic equation (11421) obtained in the 
overdamped limit becomes 

t-^, + ,/^R, W . (143) 

Introducing a mobility \x = l/(m£) and a diffusion coefficient v = D/^ 2 , this yields the stochas- 
tic equation ( l90i) without the rotational term. This can provide a microscopic model of Debye- 
Hiickel electrolytes in a canonical setting [31] which is valid in d = 2 or d = 3 dimensions. 




31 



B Correlation functions 

Considering the correlation function of the exact vorticity field fHUj) . and introducing the one 
and two-body distributions, we find that 

N 

MrH(r')) = < 7 2 J2 5 ( r - r ^ r ' " r ^> = ^ E 6 ( r " r ^ r ' ~ r )> 

i,j i=l 

+ (7 2 E 6 ( r ~ r ^( r ' ~ r J')) = ^7 2p i(r)5(r - r) + iV(iV - l)7 2 P 2 (r, r). (144) 

We now assume that the system is spatially homogeneous in average, with a vorticity distribu- 
tion uj = N^Po. Writing 5u(r) = Ud(r) — u>, we get 

(5cu(r)5cu(r')) = (u d (r)u d {r>)) - u 2 . (145) 

Starting from the identity 

(u d (r)u d (r')) = iV T 2 P 5(r - r') + N(N - l) 1 2 P 2 (r, r'), (146) 

and introducing the correlation function h(r — r') through the defining relation P2(r, r') = 
P 2 (1 + h{r - r')), we obtain for TV > 1: 

(<fo;(r)<?u;(r / )) = 7cu5(r - r') + u 2 h{r - r'). (147) 

The Fourier transform (spectrum) of the vorticity correlations is 



(27T) 2 

Finally, using the result (I37|) . we obtain 



1 + (2n) 2 -h(k) 

7 



5(k + k'). (148) 



(tektek ,> = g__L_ i(k + k0 . (149 ) 

C A heuristic relaxation equation for point vortices in 
the microcanonical ensemble 

As explained in the Introduction, the dynamical models introduced in this paper correspond 
to a canonical description of point vortices. In these models, the inverse temperature (3 is 
fixed and the energy E[t) is not conserved. Following the procedure of [441 [25| 145] , we can 
introduce heuristically a dynamical model of point vortices corresponding to a microcanonical 
description by letting the temperature f3(t) depend on time so as to rigorously conserve the 
energy. Therefore, we now assume in Eq. ([2]) that /i(t) = vf3(t) is a function of time. It is 
determined by writing E = J Pc^P/v dri...<iryv = and using Eq. ([61). This yields 

P{t) = l ~ X 8r ' 8r ' 2 — . (150) 



7 2 /Eii^v(g) d ri ...dr N 



In the mean field limit iV — > +oo, we obtain an equation of the form 



, u • Vcj = uW ■ Nuj + (3(thuVip) , (151) 
at 
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where f3(i) is determined by the condition E — J ipd t oJ dr = yielding 

^ ^ f 7u;(V^) 2 dr ^ 
We can then easily establish (see, e.g., [45]) that 

S = v I — (Vto + p-fcuWip) 2 dr. (153) 

The relaxation equation (I15ip with the constraint (11521) has the following properties: (i) T and 
E are conserved (ii) S > (iii) 5* = iff to = Ae~^ (iv) a steady state is linearly dynamically 
stable iff it is a (local) maximum of S at fixed circulation and energy. By Lyapunov's direct 
method, we know that if S is bounded from above, Eq. (I15ip will relax towards a (local) max- 
imum of S at fixed circulation and energy (if several local entropy maxima exist, the choice of 
the entropy maximum will depend on a notion of basin of attraction). Therefore, this relaxation 
equation tends to the microcanonical statistical equilibrium. This equation was introduced in 
[25] by using a maximum entropy production principle (MEPP) and it is a particular case of 
the general class of nonlinear mean field Fokker-Planck (NFP) equations introduced in [431 |4"6] . 
Note that a similar heuristic microcanonical model has been introduced for self-gravitating sys- 
tems in the form of Kramers and Smoluchowski equations with a time dependent temperature 
[TTJ [52] . We stress, however, that Eqs. ( 115ip -( fl~52~I) do not describe the real dynamics of point 
vortices in the microcanonical ensemble. The exact kinetic theory of point vortices is the one 
developed in [26] and it leads to very different equations. Therefore, Eqs. f ll51l) -( |T52l have 
probably no physical justification. However, they can be used as a numerical algorithm to 
construct statistical equilibrium states of point vortices in the microcanonical ensemble (note 
in this respect that it is possible to generalize the MEPP so as to include the conservation 
of angular momentum in a disk or in an infinite domain). Similarly, Eq. (T2"T|) can be used 
as a numerical algorithm to construct statistical equilibrium states in the canonical ensemble. 
These remarks give them some practical interest. Note also that the dynamical evolution of the 
microcanonical equation (11511) with the constraint (11521) is very different from the dynamical 
evolution of the canonical equation (T2~T|) with fixed (3. Indeed, in the microcanonical ensemble, 
there exists a statistical equilibrium state for any value of the energy E so that Eq. (11511) with 
Eq. (11521) will relax towards this equilibrium state. By contrast, in the canonical ensemble, 
there is no equilibrium state for f3 < (3 C and, in that case, the relaxation equation (f27|) leads to 
vortex collapse. 

For the multi-species system, the heuristic microcanonical model in the mean field limit is 
^ + u-Vu a = uV ■ (Vu a + (3{t) la u a V^) , (154) 
where (3(t) is determined by the condition E = J ipdtuj dr = yielding 

J cj 2 (W>) dr 

where u> = J2 a u a and u 2 = Yl a la^a- We can easily establish that 

S = vY" I ' -^—(Vuo a + (3-i a u a Vi)) 2 dr. (156) 

The relaxation equations (11541) with the constraint (11551) have the following properties: (i) 
r Q = J uj a dr and E are conserved (ii) S > (iii) S = iff uj a = A a e~ l3 ' 1a ^ (iv) a steady state is 
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linearly dynamically stable iff it is a (local) maximum of S at fixed circulations T a and energy. 
Therefore, if S is bounded from above, these equations will relax towards the micro canonical 
statistical equilibrium state of the multi-species point vortex gas. These equations can also be 
derived from the maximum entropy production principle [5T] . 
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Domain of validity 


Hamiltonian 


Brownian 


Microscopic 
iV-body model 
ror r^rj 


Deterministic Kirchhoff 
equations 
^see n/q. i^zj oi |zuj ) 


Stochastic Langevin 
equations 

(coo "Rn 

l^see n/q. i z 1 1 


Exact 
kinetic equation 
ror i^Ar^ri, r^r, rj 


Liouville equation 
(see Eq. (3) of [25]) 


iV-body 
Fokker-Planck equation 
^see ii/q. ^oj i 


Exact 
kinetic equation for 

u d (r,t) = EiIiT^( r - r i(*)) 


Klimontovich equation 
(see Eq. (69) of [26]) 


Stochastic kinetic 
equation 
(see Eq. (42)) 


Kinetic equation for 
w(r,t) = JV7Pi(r,t) 
when A 7 — * +oo 


Euler-Poisson equation 
(see Eq. (35) of [26]) 


Fokker-Planck- Poisson 
equation 
(see Eq. (25)) 


Violent relaxation 
of the 2D Euler equation. 
Kinetic equations for 
the coarse-grained 
vorticity u;(r, t) 


Phenomenological equations 
obtained from the MEPP 

(see Eq. (15) of [B] 
and Eqs. (23)-(30) of [72]) 
or kinetic equation obtained 
from a quasilinear theory 
(see Eq. (18) of [73]) 




Diffusion coefficient for 
the violent relaxation 


See Eq. (8) of [7J, Eq. (4.9) 
of [7T] or Eq. (27) of [73] 




Integra differential kinetic 
equation for ou(r,t) 
at the order 0(1/N) 


See Eq. (54) of [26] 
(collective effects neglected) 


unknown 


Kinetic equation at 
the order 0(1/ N) 
for axisymmetric flows 


See Eq. (11) of [28] 
(with collective effects) and 

Eq. (48) of [26] 
(collective effects neglected) 


unknown 


lest particle m a 
(thermal) bath 


rokker- Planck equation 
with a /ixed potential. 
See Eq. (141) of [26] 
or Eq. (140) of [2S] (thermal) 




Diffusion coefficient 
for a test particle 
in a (thermal) bath 


See Eq. (102) of [26] 
or Eq. (136) of [26] (thermal) 





Table 1: Summary of the different kinetic equations for Hamiltonian and Brownian vortices. 
A similar Table would apply to Hamiltonian and Brownian systems of material particles with 
long-range interactions [H El El [7J [8] . 
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P<Pc 


p < 


p > 


P>Pt 


2D self-gravitating 
Brownian particles 






(m, m) attract 
=>- cluster 


Collapse: 
=>- Dirac peak 
of mass Nm 


2D neutral 






(+e, +e) repel 

1 — p — p] tptipI 

(+e, — e) attract 
=>- homogeneous 
state 


Collapse: 
N/2 ^incmlar 

"atoms" (+e, — e) 


Single species 
2D Brownian 
vortices 


Collapse: 
=>• Dirac peak of 
circulation N'-f 


(7, 7) attract 
=>- cluster 
(monopole) 


(7, 7) repel 
=>- concentration 
at the boundary 




Multi species 
2D Brownian 
vortices 


Collapse: 
singular 
macroscopic 
dipole 

{+h>-h) 


(+7, +7) attract 
(—7, —7) attract 
(+7, -7) repel 
=>- two clusters 
of opposite 
sign (dipole) 


(+7, +7) repel 
(-7, -7) repel 
(+7, —7) attract 
=>- homogeneous 
state 


Collapse: 
N/2 singular 
microscopic 
dipoles (+7, —7) 



Table 2: Schematic "phase diagram" showing the analogies between 2D self-gravitating systems 
(or bacterial populations [32]), 2D Coulombian plasmas and 2D point vortices in the canonical 
ensemble [38] • For a system of N self-gravitating Brownian particles with mass m: fi+ = 
4/ '{GMm). For a neutral plasma with N/2 charges +e and N/2 charges — e: /?+ = 2/e 2 . 
For a system of iV Brownian point vortices with circulation 7: /3~ = —8ir/(N^ 2 ). For a 
multi-species system of Brownian point vortices with N/2 vortices +7 and N/2 vortices —7: 
P~ = — IQtt / (N'-f 2 ) and /3+ = 47r/7 2 . In the microcanonical ensemble, P(E) is a function of the 
energy. There is an equilibrium state for each value of the energy E, so there is no collapse. 
The critical temperatures Pf correspond to E — > =1=00 so the states ft < P~ and ft > P£ are 
not accessible in the microcanonical ensemble. 
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